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ABSTRACT 

As a pilot study for the first all-sky radio survey at short wavelengths, we have ob- 
served 1216 deg 2 of the southern sky at 18 GHz (16 mm) using a novel wide-band 
(3.4 GHz bandwidth) analogue correlator on one baseline of the Australia Telescope 
Compact Array (ATCA). We scanned a region of sky between declination —71° and 
—59° with an rms noise level of 15mJy. Follow-up radio imaging of candidate sources 
above a 4cr detection limit of 60mJy resulted in 221 confirmed detections, for which 
we have measured accurate positions and flux densities. For extragalactic sources, the 
survey is roughly 70% complete at a flux density of 126 mJy and 95% complete above 
300 mJy. Almost half the detected sources lie within a few degrees of the Galactic plane, 
but there are 123 sources with |6| > 5° which can be assumed to be extragalactic. The 
differential source counts for extragalactic sources in the range O.Uy < Sisghz < 3 Jy 
are well fitted by a relation of the form n(S) — 57 (S'/Jy) _2 2±0 2 Jy" 1 sr _1 , in good 
agreement with the 15 GHz counts published by Taylor et al. (2001) and Waldram 
et al. (2003). Over 70% of the extragalactic sources have a flat radio spectrum 
( a o 8 843 > — 0.5, S„ oc v a ), and 29% have inverted radio spectra (ao 8 8 4 3 > 0). The 
optical identification rate is high: 51% of the extragalactic sources are identified with 
stellar objects (candidate QSOs), 22% with galaxies and only 27% with faint optical 
objects or blank fields. 

Key words: surveys — cosmic microwave background — galaxies: radio continuum 
— galaxies: active — radio continuum: general 



1 INTRODUCTION 

The radio-source population above 5 GHz has not yet been 
well studied. This is mainly because large radio telescopes 
typically have fields of view of a few arcmin at high fre- 
quencies, making it extremely time-consuming to carry out 
large-area surveys. Measuring the high-frequency proper- 
ties of extragalactic radio sources is also crucial for inter- 
preting the high-sensitivity and high-resolution maps of the 
Cosmic Microwave Background (CMB) provided by WMAP 
(Bennett et al. 2003) and forthcoming missions like ESA's 
Planck satellite (Tauber 2001). Current analyses (De Zotti et 
al. 1999a) indicate that the best regions for high-resolution 
CMB mapping lie near 100 GHz, and that the main con- 
taminant is likely to be the fluctuations due to individual 
foreground radio sources. 

The brightest high-frequency radio sources are ex- 



pected to be blazars (high and low-polarization flat- 
spectrum quasars and BL Lac objects) and galaxies with 
highly-inverted radio spectra such as the GHz Peaked Spec- 
trum (GPS) sources (O'Dea 1998). The Michigan 8 GHz 
(A3. 8 cm) survey covered ~ 2600 deg 2 of sky and detected 
55 sources down to a limiting flux density of 0.6 Jy. These 
sources were found to be different from those detected in 
low-frequency surveys, with 70% having flat or inverted 
spectra. More than 80% had optical counterparts on the 
Palomar Sky Survey, most of them QSOs (Brandie & Bridle 
1974; Brandie, Bridle & Kesteven 1974). 

The first blind radio survey above 8 GHz was carried 
out by Taylor et al. (2001) with the Ryle Telescope. They 
observed 63 deg 2 of sky and detected 66 sources to a limit- 
ing flux density of 20 mJy at 15.2 GHz. A significant fraction 
(20%) of the sources they detected would not have been pre- 
dicted by a simple extrapolation from lower-frequency sur- 
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veys. Two-thirds of their sources had flat or inverted spec- 
tra. The source density observed by Taylor et al. was also 
a factor of 1.3-2 lower than predicted by theoretical models 
based on low-frequency radio luminosity functions (Dunlop 
& Peacock 1990; Toffolatti et al. 1998). 

The 15 GHz survey was extended by Waldram et al. 
(2003) to 520 deg 2 in three regions of sky targeted by the 
Very Small Array (Watson et al. 2003). They found 465 
sources above a completeness limit of ~25 mjy in their main 
survey, although a total of ~760 sources were detected down 
to a limiting flux density of 10 mjy. 

Preliminary radio source counts at 30 GHz have been 
measured by the DASI (Leitch et al. 2002), VSA (Grainge 
et al. 2003) and CBI (Mason et al. 2003) groups, and a 
catalogue of 208 point source candidates stronger than about 
1 Jy from the first year of WMAP data (five frequency bands 
from 23-94 GHz) has recently been published by Bennett et 
al. (2003). 

In this paper, we present the results of a 18 GHz Pilot 
Survey at the Australia Telescope Compact Array (ATCA) 1 . 
The pilot survey covered 1216 deg 2 of sky in the region be- 
tween 5 — —71° and —59° to a (4a) limiting flux density 
of 60mJy. In the follow-up observations we confirmed 221 
sources. 

In Section 2, we briefly describe the survey strategy 
and give some details of the observing run. Section 3 de- 
scribes the data-reduction strategy for the scan data, which 
has some novel features. Follow-up synthesis imaging obser- 
vations of the detected sources are described in Section 4, 
and their data reduction in Section 5. Section 6 presents the 
follow-up results, Section 7 a comparison with surveys at 
other frequencies, Section 8 the optical identifications and 
Section 9 summarizes our conclusions. 



2 SURVEY STRATEGY AND OBSERVATIONS 

An analogue wide-band prototype correlator originally de- 
veloped for the Taiwanese CMB instrument AMiBA (Lo et 
al. 2001) was linked to two ATCA antennas to produce a 
two-element interferometer with 3.4 GHz bandwidth. The 
two 22m antennas were placed on the shortest ATCA base- 
line (30 m) , giving a fringe spacing of ~ 90 arcsec at the cen- 
treband frequency of 18 GHz. The correlator produced 16 
lag channels corresponding (after the Fourier inversion) to 8 
frequency channels, each with a bandwidth of 425 MHz. This 
large total bandwidth allowed us to obtain high signal-to- 
noise measurements in the very short integration time used 
(0.2 sec/beam). Such rapid scanning was a design feature of 
the ATCA implemented to perform mosaic observation over 
large areas. This combination of wide bandwidth and fast 
scanning has made it possible to observe a large fraction of 
the sky despite the small ATCA primary beam at 18 GHz 
(HPBW = 2.3 arcmin). Only a single linear polarization was 
implemented in this correlator configuration. 

No delay tracking had been built for this analogue cor- 
relator (AMiBA antennas are on a common tracking mount) 
so we designed an observing mode in which all observations 

1 See www.narrabri.atnf.csiro.au and Journal of Electrical and 
Electronics Engineering, Australia, Special Issue, Vol. 12, No. 2, 
June 1992 



are taken on the meridian where the delay for an east-west 
interferometer is zero. 

The scanning strategy was as follows: the two ATCA 
antennas point in the same direction and scan the sky syn- 
chronously up and down between 8 — —70° and —60° in ele- 
vation along the meridian with a scan rate of 10° per minute 
(the maximum speed the antennas can slew in a closed loop 
in elevation). The Earth's rotation then allows us to cover 
a strip of the sky in raster scan mode. Antenna coordinates 
were sampled at 10 Hz. With an integration time of 80 ms, 
the antenna pointing was recorded every 0.8 arcmin. 

The pilot survey run took place on 2002 September 13- 
17 during one observing period of 96 hours. 



3 SURVEY DATA REDUCTION 

The antenna positions (obtained from the telescope en- 
coders) and fringe visibilities (produced by the analogue cor- 
relator) were stored in data files. Custom MIRIAD (Sault, 
Teuben & Wright 1995) tasks were written to flag out bad 
data (e.g. antenna pointings more than 30" from one an- 
other) and to calibrate the visibilities in absolute flux and 
bandpass. 

The next stage involved searching for peaks in the fringe 
visibility amplitudes above a given signal-to-noise ratio in 
the time-ordered raster scans, to produce a candidate source 
list. This was done using the following steps: 

(i) We calculated the weighted complex mean of the 8 sep- 
arate frequency channels to form a single 3.4 GHz wide band. 
The bandwidth smearing (Thompson 1999) in the worst case 
(a source at half power point of the primary beam) is 1/6 of 
the fringe spacing and can be ignored. Expressed in a differ- 
ent way our 15% bandwidth corresponds to a 5.6m variation 
in the 30m spacing used, which is a small fraction of the 22m 
dish diameter. 

(ii) We apply a beam shape Gaussian (FWHM = 2'. 3) 
filter acting as a weighted running mean for this short base- 
line to smooth out fluctuations (generated by receiver noise) 
on spatial scales less than the antenna primary beam. 

(iii) In the plane defined by the real an imaginary parts 
of the fringe visibility we dete rmined t he dispersion along 
each axis. Points lying at > n^J <r 2 + Oy are selected as can- 
didate sources at na level. In general a source generates up 
to 4 contiguous points along the time series corresponding 
to the sampling by different parts of the beam in steps of 
0.8 arcmin. The candidate source position corresponds to the 
maximum visibility amplitude in this sequence. 

(iv) We created a list of candidate point-source detec- 
tions at a range of signal-to-noise ratios (3, 3.5, 4 and 5 a) 
to explore the optimum detection strategy. 

The raw candidate source lists were then merged into 
a single list. The number of 4a candidate sources for the 
entire survey was 574. The measured rms noise in the scans 
was ffrms = 15mJy. 

We performed a Monte Carlo simulation to estimate the 
number of spurious peaks generated by Gaussian noise at a 
given signal-to-noise ratio (see Table 0. The same number 
of real and imaginary parts of the fringe visibilities for each 
frequency channel was generated 30 times by sampling a 
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N cxp 


N aim 


3 


542 


529±21 


3.5 


in 


104±10 


4 


24 


16±4 


5 


5 


0.13±0.35 



Table 1. Comparison between the actual number of peaks in a 
raster scan during one test run (N oxp ), and the number of spurious 
peaks predicted by a Monte Carlo simulation of the same scan 
(N s ; m ) as described in the text. 



Gaussian distribution with the same mean and standard de- 
viation as the actual distributions. A cutoff value of S/N = 4 
was chosen as a good trade-off, since lower values yield too 
many spurious peaks while for S/N = 5 we reject too many 
sources which are genuine detections. Since the S/N in the 
follow-up is very high compared with the 0.2 sec survey in- 
tegration, we could afford to explore the detection statistics 
by going well into the region of spurious detections due to 
noise. 



4 FOLLOW-UP OBSERVATIONS 

Follow-up synthesis-imaging observations of the S/N > 4 
candidate sources were carried out on 2002 October 8-12 
with the ATCA in its hybrid H168B configuration, which has 
antennas along the new 214 m Northern Spur as well as the 
existing east-west railway track. At the time of this run, only 
three antennas had been equipped with low-noise 12 mm re- 
ceivers, so only three baselines were available for imaging. 
We therefore observed each object over a wide range in hour 
angle to maximize uv coverage. The angular resolution de- 
pends on the source declination. The median major and mi- 
nor axis values are 30 and 13 arcsec respectively. The average 
rms noise in synthetized maps was a ~ 1.5 mjy. 

We used the standard ATCA correlator, allowing us 
to observe in two spectral windows centred at 17.2 and 
18.7 GHz each with a bandwidth of 128 MHz. At 18 GHz, the 
standard ATCA primary calibrator 1934 — 638 is too faint to 
provide a good absolute flux density scale, and we used the 
planet Mars as our primary calibrator. The fringe-visibility 
amplitude dependence on spatial frequency has been accu- 
rately modelled for the ATCA. The Mars model we used 
is described in Ulich (1981) and Rudy (1987). As in Ulich 
(1981), the seasonal variation of Mars brightness tempera- 
ture is also taken in account by the model. 



5 FOLLOW-UP DATA REDUCTION 
5.1 Follow— up 18 GHz images 

Reduction of the follow-up 18 GHz images was done in the 
standard way using the MIRIAD software package (Sault, 
Teuben & Wright 1995). The target sources and calibra- 
tors were imaged in total intensity (Stokes I) and circular 
polarization (Stokes V). The sources were assumed to have 
zero circular polarization, and the Stokes V maps were used 
to estimate the rms noise. We used the MIRIAD Multi- 
Frequency Synthesis (MFS) algorithm at the Fourier inver- 



sion stage, since this improves (u, v) coverage by simulta- 
neously Fourier-transforming the visibility data from both 
observing frequencies. After the Fourier inversion the maps 
were CLEANed, and the resulting clean components were 
restored using an elliptical Gaussian beam. 

We used the maximum pixel value in the restored Stokes 

1 map together with the V Stokes rms noise to calculate 
the map signal-to-noise (hereafter S/N) ratios. A source was 
considered to be detected when its S/N value was above 5. 

The number of sources above this detection threshold 
was 282 out of 574 candidates. However, many candidate ob- 
jects (~ 20%) were multiple detections of the same source, 
so that the actual number of spurious sources in the sur- 
vey was smaller. These multiple observations were used to 
estimate data quality. 

5.2 Flux density and position measurements 

We reimaged the detected sources in all four Stokes parame- 
ters, applying the primary beam correction to restored maps 
and a natural weighting to the fringe visibilities before the 
Fourier inversion. The primary beam correction was only 
applied at this last stage, because the increase in noise near 
the map border complicates the detection statistics. 

All but five of the confirmed (S/N> 5) sources were 
well-fitted by a single elliptical Gaussian. 

We created a final source list by removing all sources 
which had the same position (within the errors) as a 
previously-listed one. This decreased the final number of 
detected sources from 282 to 221, and most of the subse- 
quent analysis was done on this list of objects. Source pa- 
rameters determined by the MIRIAD IMFIT algorithm 2 for 
the confirmed sources are presented in Table 2. Any source 
having an angular size larger than 10 arcsec was dubbed 
(R) meaning "significantly resolved" . Because of angular size 
measurement errors, point-like sources at the resolution of 
follow-up images could appear resolved. 

5.3 Polarization 

Determining the level of polarization in foreground radio 
sources is critically important for the study of polarization 
anisotropies in the cosmic microwave background (e.g. Hu, 
Hedman & Zaldarriaga 2003). De Zotti et al. (1999b) and 
Mesa et al. (2002) estimate that extragalactic sources will 
not be a strong limiting factor in measurements of E-mode 
CMB polarization, but there are few existing measurements 
of the polarization of radio galaxies and quasars at frequen- 
cies above 5 GHz (Ricci et al. 2004). 

For all the sources and calibrators we produced polar- 
ized intensity p, fractional polarization m and position an- 
gle PA maps and the relative error maps after clipping (i.e. 
blanking) pixels with S/N ratios below 3. The peak flux in 
the Q, U and I maps was then measured at the I map peak 
position using the MIRIAD MAXFIT task (see the MIRIAD 
web page in the previous footnote) , and the fractional linear 
polarization estimated through the relation: 

2 Described in wwwatnf.csiro.au/computing/software/miriad/ 
taskindex.html 
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Table 3. Linear polarization detections at 3 a. 



v/Q 2 + U 2 



(1) 



Only 13 of the 221 follow-up confirmed sour ces were 
found to have a linear polarization intensity p = */ Q 2 + U 2 
above 3a p , where a p is typically 6mJy (Table [3J. Six of 
these have |fc| > 5° and appear to be extragalactic. Two 
extra-galactic objects in the \b\ < 5° Galactic strip were 
also detected in polarization. 

In Fig. we show the fraction of 18 GHz polarized 
sources that were actually detected below a certain polariza- 
tion level inside the sample of sources that could be detected 
below that polarization level. This statistical method, simi- 
lar to the one used by Auriemma et. al (1977) to study ra- 
dio luminosity functions, gives an estimate of the fractional 
polarization distribution. Error bars indicate the maximum 
range inside each polarization sub-sample. The first polar- 
ization level is strongly affected by small number statistical 
Poissonian errors, having 7 polarization detections out of 8 
sources in which polarization could be detected. If we take 
only the extragalactic sources into account we have a me- 
dian polarization of ~ 2.1%. The average flux density of all 
extragalactic sources contributing to the median is 1.0 Jy. 

We also draw attention to the extragalactic source 
1824-6717 (=PKS 1819-67), which has the highest frac- 
tional polarization so far measured in our sample at 16.2%. 
Since PKS 1819— 67 is a widely-spaced (47arcsec) double, 
with extended radio emission seen at lower frequencies 
(Fig. |2l , it is unlikely to be a distant gravitationally-lensed 
radio galaxy. We are probably seeing a pair of radio hotspots 
in a powerful FRII radio galaxy. Since such sources are iden- 
tified with intrinsically bright elliptical galaxies and no op- 
tical counterpart is seen on the DSS, the galaxy redshift 
should be z > 1. 



6 DISTRIBUTION OF SOURCES 

6.1 Distribution of the detected sources on the 
sky 

As can be seen from Figs. [3] and 2] many of the confirmed 
18 GHz sources lie in or near the Galactic plane. 

Using the source density histogram in Fig. |3] we esti- 
mated the number of galactic objects that are expected to 
be found in the extragalactic survey region. The galactic 




Figure 1. Fractional polarization distribution at 18 GHz given 
by the ratio between the number of polarization detected sources 
and the number of sources that could be detected at this level. 
Errors bars provide the maximum range permitted by Poissonian 
errors in source sample size per bin. 



PKS 1819-67 




Figure 2. Radio contours for PKS 1819-67 at 843 MHz (larger 
scale) and 18 GHz (smaller scale) overlaid on a blue Supcr- 
COSMOS image. This source has the highest linear polarization 
(16.2%) measured in our pilot study. 

plane shape in the source density histogram was fitted with 
a single Gaussian profile. The best fit parameters resulted 
to be n pe ak=2.3 deg -2 , &o=0 and cr=0°.8, where n peak is the 
source density peak, bo is the galactic latitude value at the 
source density peak and a is the Gaussian profile width. 

In the discussion which follows we assume that the 
123 sources which have |6| > 5°, and which lie more than 
5.5 degrees from the LMC centre (a = 5 h 23 m 34 s .7, 6 = 
-69° 45' 22", J2000; I = 280.47°, b = -32.89°) are extra- 
galactic. When split in this way, the final sample of con- 
firmed sources in Table 3 contains 84 Galactic, 14 LMC and 
123 extragalactic sources. 
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Figure 3. Histogram of source density in 1-degree bins of Galac- 
tic latitude b (corrected for the area surveyed). The peak near 
6=0° corresponds to Galactic disk sources such as HII regions. 
The surveyed area (see Fig. |3J only extends to b = +3° so the 
apparent cut-off in sources at b > +3° is not real. 
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Figure 5. Effective area A er T surveyed in the extragalactic region 
of the pilot survey scans as a function of the limiting flux density 
(•^lim = 4 X AS rms ). The maximum surveyed area outside the 
Galactic latitude range (|6 | > 5°), excluding the LMC region, is 
A tot = 1130 deg 2 . 




Figure 4. Equal area Aitoff projection sky map in galactic co- 
ordinates of the confirmed 18 GHz sources. Larger crosses corre- 
spond to brighter sources. 



6.2 Survey sensitivity 

Because this pilot survey did not achieve uniform cover- 
age we made a quantitative assessment of the survey sky 
coverage through a survey sensitivity map, i.e. a sky map 
representing the limiting flux density value at every point 
of the sky covered by the survey. To create this map, the 
survey sky strip (—71° < S < —59°) was gridded into 
30 arcsecx30 arcsec pixels, and the ~ 2736890 survey sample 
pointings were projected onto this 1800x24000- pixel grid. 
To account for the antenna primary-beam pattern, this im- 
age was convolved with a 7 x 7-pixel discretized Gaussian 
function (HPBW = 2.3arcmin). The convolved pointing 



number per pixel was then multiplied by the time t CTOSS taken 
by a single beam to sweep across a pixel (icross = 0.28 sec) 
to get the pixel observing time. The antenna system tem- 
perature (T sys = 78 K), the antenna gain at 12mm (G = 
0.10 K Jy _1 )i the observing time per pixel r and the corre- 
lator band-width (Av — 3.4 GHz) were used to compute the 
pixel rms sensitivity ASVins- A quantitative representation 
of the sensitivity map is offered by the effective area A e g 
as a function of the limiting flux density (Fig. |SJ, which is 
simply given by the cumulative number of pixels for which 
4 x ASrms < Sum multiplied by the single pixel area. 



6.3 Differential counts of extragalactic sources 

Here, we determine source counts for the extragalactic pop- 
ulation at 18 GHz. We exclude Galactic sources, since many 
have large angular size and their integrated flux densities 
are underestimated in the follow-up snapshot images. 

The counts were logarithmically binned in flux density 
starting from Sn m =100mJy with a common binsize of 0.2 
dex. A limiting flux density of 100 mjy was chosen because 
only a handful of objects were detected below this flux den- 
sity. 

As the effective surveyed area A c ff drops with decreas- 
ing limiting flux density Sii m below Sum — 1 Jy, the sur- 
face density of sources in the (small) flux density interval 
AS* = Sa — Si, n(S)AS, is given by: 



n(S)AS = N obs (> Si) 



A cS (> Si) 



-N obs (> Sa) 



1 



A eH (> Sa) 



,(2) 



where N h s (> S) is the number of counted sources brighter 
than S. This procedure is equivalent to the more commonly 
used method of calculating the effective area over which each 
source could have been observed and summing the inverse 
area in the flux density bin of interest (Katgert et al. 1973). 
The computed source counts need to be corrected for 
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Table 5. Differential counts of extragalactic sources, corrected 
for incompleteness. 
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incompleteness. The fraction of sources with true flux S > 
4a, that should be in our sample but are missed because 
their observed flux Sobs < 4<j is given by: 



S/4 pitr 

da A(a) I dS ohB e 
™„ Jo 



2 V " I /V2^a 



S/4 



.(3) 



da A(a) 



The values of Flnissed as a function of S are shown in TableQ3 
together with their negative and positive errors. 

We have used the lowest possible angular resolution 
in the survey to minimize the selection against extended 
sources, but any sources extended more than a few minutes 
of arc could be missed at the survey limit. This will have a 
very small effect on the estimate of the total source counts. 

The uncertainties on differential counts were estimated 
by summing in quadrature the Poisson errors (Gehrels 1986) 
and the uncertainties in the effective area. The results are 
plotted in Fig. |S| where we also show the least square linear 
fit given by: 



Jy sr~ 



(4) 



Our counts are in good agreement with those by Waldram 
et al. (2003) (9C survey) at 15 GHz (Fig.Q, as well as with 
the preliminary VSA (Grainge et al. 2003), DASI (Leitch 
et al. 2002), and WMAP counts (Bennett et al. 2003) if we 
assume a radio spectral index a = 0. 
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Figure 6. Effective area-corrected and completeness-corrected 
differential source counts of the extragalactic population at 18 
GHz are presented together with 9C main survey (triangles) and 
deeper survey (squares) counts by Waldram et al. (2003). The 
solid line shows the least-squares linear fit to our counts. The 
two lowest flux density bins (triangles) of the 9C main survey are 
underestimated due to incompleteness. 
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7 COMPARISON WITH DATA AT LOWER 
RADIO FREQUENCIES 

7.1 Position comparison 

Twenty of the objects detected in the 18 GHz pilot survey 
are strong compact radio sources which have previously been 
observed by Ma et al. (1998) with Mark III VLBI at 2.3 and 
8.4 GHz as part of a program to tie down the International 
Celestial Reference Frame. The mean offset is +0.5arcsec in 
RA and —0.2 arcsec in Dec, with standard deviations of 3.0 
and 1.4 arcsec respectively. These standard deviations are 
larger than the formal uncertainties (2.0 arcsec in RA and 
1.1 arcsec in Dec) measured in the Miriad source-fitting, and 
reflect the true external errors in our position measurements. 

We also cross-matched our extragalactic sub-sample 
with the 0.843 GHz SUMSS (Mauch et al. 2003) and 



Figure 7. Euclidean normalized differential source counts from 
our pilot survey (18PS) compared with VSA, DASI, WMAP band 
Q (41 GHz) and 9C main (triangles) and deeper (squares) surveys. 
Modified version of Fig. 13 from Bennett et al. (2003). 



4.85 GHz PMN (Gregory et al. 1994) catalogues. All the 123 
18 GHz extragalactic sources have a PMN counterpart, and 
all 102 sources lying in the region of overlap with SUMSS 
have a SUMSS counterpart. 

7.2 Spectral index distribution 

Fig. |H| shows that the low-frequency (So. 843) and high- 
frequency (Sis) flux densities are uncorrelated. It is im- 
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Figure 8. So. 843 vs. Sis scatter plot. No evidence of a significant 
correlation between low- and high-frequency fluxes is present. 
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Figure 9. Spectral index distributions of the extragalactic sam- 
ple. The solid line shows a histogram of the spectral indices ajj g43 
calculated from SUMSS at 843 MHz (median value=-0.22), and 
the dashed line «g 8 calculated from PMN at 5 GHz (median 
value=-0.29). 



portant to stress this point because in many CMB exper- 
iments the extra-galactic radio source foreground is com- 
puted starting from low-frequency source catalogues and 
assuming a distribution of spectral indices around a mean 
value to model the source population at high frequency (>30 
GHz), but in doing this, a correlation between low- and high- 
frequency flux density is implicitly assumed. 

Fig. |3] shows the spectral index distribution (a where 
S v oc v a ) for the extragalactic population. Most spectral 
indices are in the range —1.5 < a < 0.5. Over 70 % of the 
extragalactic sources have a flat radio spectrum (ao 8 843 > 
—0.5) and 29% have an inverted spectral index (ofJ 8 843 > 0). 

The cvo 8 843 spectral index distribution we obtained is 



in good agreement with the Michigan 8 GHz survey (70% 
flat or inverted sources) and Taylor et al. (2001) 's distribu- 
tions (2/3 of the sources flat or inverted). However, the 9C 
(Waldram et al. 2003) spectral index distribution between 
1.4 GHz NVSS and 15 GHz flux densities, which goes about 
4 times fainter, contains almost equal numbers of fiat and 
steep spectrum sources (55% flat or inverted), suggesting 
that the ratio of flat to steep spectrum sources decreases at 
lower flux density. The same behaviour is found for the spec- 
tral index statistics of a lower frequency survey such as the 5 
GHz selected 1 Jy sample (54% of flat or inverted spectrum 
sources, Stickel et al. 1994). Theoretical predictions based 
on low frequency selected radio luminosity functions (e.g. 
Dunlop & Peacock 1990) and simple spectral scaling laws 
(flat spectrum a = and steep spectrum a = —0.7) tend 
to overestimate the source density at high frequency and 
fail to reproduce the observed proportionality between flat 
and steep spectrum populations. High-frequency turn-over 
in radio source spectra due to synchrotron electron ageing 
needs to be modelled as we are at reasonably high observing 
frequency and some sources lie at high redshift. 



To better characterise the spectral behaviour of the 
extra-galactic population we also present an al, 8 vs. ag.g 
scatter plot in Fig. llOl for the 101 extragalactic objects hav- 
ing a flux density measurement at 0.8, 5 and 18 GHz. 



We found that 15 sources classified as inverted at low 
frequency (0 < Qg.g < 0.5) are also inverted at high fre- 
quency (CV5 8 > 0). No sources out of the 6 classified as High 
Frequency Peaker candidates having «o 8 > 0.5 show a rising 
spectrum between 5 and 18 GHz. 14 sources being flat be- 
tween 0.8 and 5 GHz (—0.5 < ag.g < 0) become inverted at 
higher frequency. Most of the steep spectrum («q 8 < —0.5) 
sources (except 2 objects) stay steep even in the high fre- 
quency range. 

In the «5 8 histogram in Fig.|U]there are 9 sources in the 
0.4-0.6 bin but only 2 are present in the same spectral range 
in the ajj.8 histogram. The radio spectra of these objects are 
turning up at higher frequency. It is worth noticing that 
the PMN flux densities could be overestimated as the PMN 
survey beam was so large (3 arcmin) that it is possible that 
more than one source is included inside the beam area. This 
would flatten the 0.843-5 GHz spectral indices and steepen 
the 5-18 GHz ones. 

Within the boundaries in Fig. 0] there are 12 WMAP 
sources. We found that six of them fall in small survey 
"holes" where the survey sensitivity was too poor or array 
antennas were unable to sample sky positions and ampli- 
tudes. The effect of these regions with poor or null sensitiv- 
ity on the survey effective area and completeness is quantita- 
tively taken in account in Fig. [S] For the six sources having 
a S/N > 4 in the survey we were able to produce wide range 
(0.8-61 GHz) radio spectra which in four cases turned out 
to be steeply rising spectra. One source (PMN J 1803-6507) 
could be a HFP candidate possibly peaking at about 40 
GHz (Fig.inj. Note that WMAP and 18 GHz pilot survey 
flux densities were measured within a few months of each 
other (i.e. they are near- simultaneous) , the PMN (5 GHz) 
and SUMSS (0.843 GHz) measurements date from 10 and 4 
years earlier respectively. 
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Figure 10. ctj 8 vs. g scatter plot: the dashed line represents 
the separation between flat (oq g > —0.5) and steep spectrum 
sources and while the dotted line is the threshold for selecting 
High Frequency Peaker (a^ g > 0.5) candidates at low frequency. 
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Figure 11. Radio spectrum spanning between 0.843 and 61 GHz 
of one of the 6 sources in the 18 GHz pilot survey present also 
into the WMAP discrete source catalogue. The WMAP points 
are represented by diamonds. 



of this search for the 123 extragalactic sources in our sam- 
ple. Of these, 74 (60%) are identified with stellar objects 
(candidate QSOs), 30 (25%) with galaxies and 19 (15%) are 
blank fields, i.e. they have no optical counterpart in Super- 
COSMOS. Some notes on individual sources are given in 
Appendix A. 

The table columns are as follows: 
(1) 18 GHz source name, as in Table 2. 

(2,3) 18 GHz radio position (J2000.0) measured from the 
follow-up images as described in §5.1. 

(4,5) The position (J2000.0) of the closest optical coun- 
terpart identified in the SuperCOSMOS catalogue. Where 
these columns are blank, there is no catalogued optical 
object within 8arcsec of the 18 GHz radio position. 

(6) Blue (Bj) apparent magnitude of the optical ID. 

(7) SuperCOSMOS T parameter, where T=l corresponds 
to a spatially-resolved galaxy and T=2 to a stellar object. 

(8) Offset (in arcsec) between the radio and optical posi- 
tions. 

(9,10) Total (integrated) 18 GHz flux density and estimated 
error. Note that these are generally higher than the peak 
flux densities listed in Table 2. 

(11,12) 5 GHz flux density (in mjy) from the PMN cata- 
logue (Gregory 1994), and estimated error. 
(13,14) 843 MHz flux density (in mjy) from the SUMSS 
catalogue (Mauch et al. 2003), and estimated error. Where 
no flux density is listed, the source lies in an area of sky 
where SUMSS observations had not been completed by 
mid-2003, and which was therefore not included in the 
Jul-1-2003 version of the SUMSS catalogue. 
(15,16) Radio spectral index between 843 MHz and 18 GHz, 
and estimated error. 

(17) Source name in the NASA Extragalactic Database 
(NED). 

(18) Notes on NED sources. 

Sixteen of the objects in Table 4 have published red- 
shifts, and these range from z=0.014 to 0.183 for the galaxies 
and .z=0.54 to 3.15 for the QSOs. 

Fig. H2l shows the relation between 18 GHz flux density 
and spectral index for sources identified with QSOs, galaxies 
and blank fields. There appears to be a clear similarity be- 
tween QSO and galaxy spectral index distributions, while in 
low-frequency surveys these populations have a very differ- 
ent behaviour. We suggest that our sample is dominated by 
AGN in galaxy nuclei rather than radio galaxy lobes. This 
then points to a close similarity between radio loud QSO 
and galaxy nuclei populations. 



8 OPTICAL IDENTIFICATIONS OF 
EXTRAGALACTIC SOURCES 

8.1 Optical identifications from SuperCosmos 

We used the SuperCOSMOS online database (Hambly et 
al. 2001) to search for optical counterparts to the 18 GHz 
sources detected in the pilot study. Our radio position errors 
are typically 3.0 arcsec in RA and 1.4 arcsec in Dec (§7.1). 
The optical positions listed in the SuperCOSMOS catalogue 
also have an associated uncertainty of up to 1 arcsec, so we 
chose a cutoff radius of 8 arcsec around the radio position 
to search for an optical counterpart. Table 4 lists the results 



8.2 Reliability of the optical IDs 

We used a series of Monte Carlo tests to estimate the num- 
ber of spurious SuperCOSMOS radio-source IDs which arise 
from the chance superposition of a radio source and a fore- 
ground star or galaxy. To do this, we offset the position 
of each radio source by an amount (typically 10-30 arcmin) 
which is significantly larger than our 8 arcsec matching ra- 
dius, but small enough that the surface density of Super- 
COSMOS stars and galaxies should be similar. We then 
cross-matched the offset radio positions with the SuperCOS- 
MOS catalogue in the same way as was done for the genuine 
radio positions. This was done four times, with offsets to 
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Figure 12. Plot of 0.843-18 GHz spectral index versus 18 GHz 
flux density for sources identified as QSOs (filled circles) and 
galaxies (open circles), and for sources which are blank fields 
(Bj > 22 mag, triangles). The error bar in the bottom right hand 
corner is representative of the typical uncertainty in the spectral 
indices which is dominated by the typically 20% errors of 18GHz 
flux densities. 



the north, south, east and west of the original positions, to 
measure the average number of spurious matches. 

We found that for a cutoff radius of 8arcsec and an 
optical magnitude Bj < 22 mag, roughly 17% of the optical 
IDs in Table 4 are likely to be spurious, giving an overall 
reliability of ~83%. This could be improved by measuring 
more accurate radio positions for the 18 GHz sources, so that 
the cutoff radius for matching could be reduced. The small 
number of very faint optical IDs (Bj = 22 — 23 mag) in Table 
4 is close to the number expected by chance, suggesting that 
almost all of these very faint objects are chance associations 
rather than genuine IDs. 



8.3 Optical identification rate 

Restricting the sample to a limiting magnitude of Bj = 22.0 
where the SuperCOSMOS catalogue is complete and the re- 
liability of the radio-optical matching high, we find an over- 
all optical ID rate of 80±8% (98/122 extragalactic sources). 
When corrected statistically for chance coincidences with 
foreground objects, this drops to 73% with an optical ID 
brighter than Bj = 22.0, where 51% are stellar objects (can- 
didate QSOs), 22% galaxies and 27% faint objects or blank 
fields. This rate is comparable with the 83% ID rate reported 
by Snellen et al. (2002) for a sample of flat -spectrum ra- 
dio sources selected at 6 cm and with flux densities above 
200 mJy. It is however much higher than that found for 
weak low-frequency samples (Windhorst, Kron & Koo 1984; 
McMahon et al. 2002). 



9 CONCLUSIONS 

The wide-band interferometer (total bandwidth spanning 
from 16 to 20 GHz) produced by linking a wide-band cor- 
relator to two ATCA antennas equipped with low-noise re- 
ceivers allowed us to survey the sky at a rate of more than 
12deg 2 per hour to a detection limit of ~ 100 mJy despite 
the small primary beam area (FWHM = 2.3 arcmin). A total 
of 221 sources have been confirmed by follow-up observations 
and had flux densities and positions measured. 

The distribution with Galactic latitude shows a sharp 
peak with FWHM in |6| of 1.9° due to Galactic sources. 

We exclude the galactic plane and use only the 123 
sources with \b\ > 5°. Although this pilot survey is not ex- 
tremely uniform we know the sensitivity map and can es- 
timate the 18 GHz counts of extragalactic sources down to 
~ 100 mJy. The results are in good agreement with the 15 
GHz counts by Taylor et al. (2001), with the preliminary 
VSA (Grainge et al. 2003) and DASI (Leitch et al. 2002) 
counts, with the WMAP counts at 41 GHz (Bennett et al. 
2003) and with the 9C survey at 15 GHz by Waldram et al. 
(2003). We determine the source counts particularly well in 
the crucial 200 - 800mJy flux range, bridging the gap be- 
tween the 9C and the WMAP surveys. With our estimate of 
the slope as S -2 ' 2 we can confirm that the Planck mission 
will only need to measure and correct for foreground sources 
down to the lOOmJy level in order to push the fluctuations 
due to extragalactic radio sources below the instrumental 
noise (cf. Fig. 6 of De Zotti et al. 1999a). At this flux level 
a blind survey of the entire sky at 20GHz is feasible. 

Although this survey has not revealed any new class 
of radio sources, we drew attention to the huge scatter in 
the flux - flux scatter plot (Fig. |HJ which suggests that care 
is needed when using low-frequency surveys to correct high 
frequency CMB observations. 

The radio spectral index distribution of the extra- 
galactic sources shows that over 70% of them are "flat"- 
spectrum (aj 8 843 > -0.5, S u oc v a ), including 30 (29%) 
sources with rising spectrum (ao 8 g43 > 0) . We have a signif- 
icant population of sources with minima in their spectrum 
in the 1-10 GHz range. 

The optical identification rate of 73% is very high: 51% 
of the extra-galactic sample are identified with stellar ob- 
jects (candidate QSOs) and 22% with radio galaxies. This is 
almost a factor of 3 higher than the identification rate in sur- 
veys at comparable flux levels at 1.4GHz. This indicates that 
the high-frequency source population has a luminosity func- 
tion which does not extend to such high power as in lower 
frequency samples. We thus expect that high-frequency sur- 
veys will be more sensitive to clustering effects. 

A statistical estimate of fractional polarization distribu- 
tion yielded a median value of 2.1% which is low compared 
with the steeper spectrum population. It agrees with the 
median polarization measured for the flat-spectrum popula- 
tion in the Kiihr sample (Ricci et al. 2004) again indicating 
that this high-frequency survey has not revealed any new 
population of sources at this flux level. With the median 
polarization below that expected for CMB E-mode, Planck 
will have no difficulty correcting for this foreground popu- 
lation. However B-mode CMB polarization is much harder 
to detect and correction for polarised sources will be neces- 
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sary. Future targeted follow-up with much higher polariza- 
tion sensitivity is possible. 

In the near future we plan to cover 10 4 deg 2 in the 
southern sky by using an improved and more powerful ver- 
sion of the wide band correlator (3 baselines and 8 GHz 
bandwidth) being able to reach a 100% completeness for all 
sources with S20GHZ > 40 mjy and we expect to detect ~ 10 4 
extragalactic sources down to that limiting flux density. 

This survey has confirmed a lower than expected den- 
sity of sources at high frequency making the mm array 
calibration problem more difficult, however the full 18GHz 
ATCA survey will still be able to provide a superb calibra- 
tion source network for ALMA with a density of about one 
source per square degree down to 40mJy. 
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APPENDIX A: NOTES ON INDIVIDUAL 
SOURCES 

PKS 0021-686: We adopt the SuperCOSMOS classifica- 
tion as a stellar object, but note that the 2MASS Extended 
Objects catalogue (Jarrett et al. 2002) classified this as 
a galaxy with K s (total) = 13.65 mag. Detected as an X-ray 
source in the RASS Bright Source Catalogue (Voges et al. 
1999). 

PKS 0022— 60: Compact double source shown in Fig. 
Al. The catalogued MRC flux density of 3.83 (±0.17) Jy at 
408 MHz implies a radio spectral index a?g 408 = -0.78±0.02. 

PKS 0101-649: The 843 MHz flux density is from 
Jones & McAdam (1992). Listed in the 2MASS Extended 
Objects catalogue as a galaxy with K s (total) = 13.41 mag. 

PKS 0235—618: Detected as an X-ray source in the 
RASS Bright Source Catalogue. Jones & McAdam (1992) 
note that the original Parkes detection is resolved into a 
pair of unrelated sources at 843 MHz. The source detected 
at 18 GHz is associated with the western (brighter) member 
of the pair. 

PMNJ 0422-6507: The optical ID is offset 7.9 arcsec 
in RA from the 18 GHz radio position, but is within 
0.5 arcsec of the SUMSS position. We accept this as the cor- 
rect ID because of the close agreement with the SUMSS 
position, and because the 18 GHz source is close to the edge 
of the observed field. 

PKS 0516-621: The optical ID is within 0.1 arcsec of 
the VLBI position measured by Ma et al. (1998). 

PKS 0522-611: Noted as a ROSAT X-ray source by 
Brinkmann, Siebert & Boiler (1994). 

PKS 1105-680: The optical ID found in SuperCOS- 
MOS was first identified by Jauncey et al. (1989), and is 
within 0.1 arcsec of the Ma et al. (1998) VLBI position. The 
foreground Galactic extinction in this field (b = —7.4°) is 
A B = 0.99 mag. 

PKS 1133-681: There is no optical object within 
5 arcsec of the 18 GHz radio position in the SuperCOSMOS 
catalogue. White et al. (1987) identify this source with a 
faint (B=22.2mag.) galaxy which is just visible on the Su- 
perCOSMOS B image. The foreground Galactic extinction 
in this field (b = -6.6°) is A B = 2.0 mag. 

PKS 1420-679: The SuperCOSMOS image shows no 
optical object within 5 arcsec of the radio position. White 
et al. (1987) identify this source with a faint (B=22.2mag.) 
stellar object, but this should be regarded as uncertain with- 
out spectroscopic confirmation because of the high surface 
density of foreground stars at this low Galactic latitude. The 
Galactic extinction in this field (6 = —6.8°) is Ab = 2.4 mag. 

PKS 1448-648: The SuperCOSMOS image shows no 
optical object within 5 arcsec of the radio position. White 
et al. (1987) identify this source with a faint (B=22.0mag.) 
galaxy which is 7.5 arcsec from the 18 GHz position, but this 
object is classified as stellar by SuperCOSMOS and may be 
a foreground star. The extinction in this low-latitude region 
(b = -5.1°) is A B = 3.2 mag. 

WKK 5585: Detected as an X-ray source in the RASS 
Bright Source Catalogue (Voges et al. 1999). 

NGC 6328: This nearby (z=0.014) galaxy is the opti- 
cal counterpart of the radio galaxy PKS 1718—649, which is 
the closest-known GPS radio source (Tingay et al. 1997). 
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PKS 1801-702: A bright elliptical galaxy without a 
catalogued redshift, K s (total) = 13.65 mag in the 2MASS Ex- 
tended Objects catalogue. 

PKS 1814-63: A nearby (z=0.0627) active galaxy, 
with K s (total) = 13.65mag in the 2MASS Extended Objects 
catalogue. 

PKS 1819—67: Highly-polarized double source with no 
optical ID, discussed in §5.3. 

PKS 1922-62: This source is 4.5 arcmin from the 
strong radio source PKS 1922—62. Large pointing errors in 
PMN (Parkes telescope beam width is about 5 arcmin) are 
likely to be the cause of the offset. 

PKS 1934-638: This strong and well-studied radio 
galaxy was the prototype GPS source (O'Dea et al. 1991). 

PMNJ 2008-6110: Compact source with a highly- 
inverted radio spectrum and no optical ID. 

MRC 2041—617: Detected as an X-ray source in the 
RASS Bright Source Catalogue (Voges et al. 1999) and 
shown in Fig. Al. 

ESO075— G41: The two components detected here 
(Fig. Al) are the core and southern hotspot of this nearby 
(2=0.028) radio galaxy. 

PKS 2210—637: Detected as an X-ray source in the 
RASS Bright Source Catalogue (Voges et al. 1999). 

PKS 2300-683: Noted as a ROSAT X-ray source by 
Brinkmann, Siebert & Boiler (1994). 

PKS 2303—656: There is no optical counterpart in the 
SuperCOSMOS catalogue, but Jackson et al. (2002) report 
that a faint red galaxy is visible on CCD images of this 
field, with a tentative redshift of 2=0.470 from an optical 
spectrum. 

IRAS 23074-5957: Noted by Roy & Norris (1997) as 
a member of the rare class of 'radio-excess infrared galax- 
ies'. Listed in the 2MASS Extended Objects catalogue with 
K s (total) = 14.00mag. 

PKS 2356-61: Well known FRII radio galaxy (Koeke- 
moer & Bicknell 1998)- both hotspots are detected at 
18 GHz. The host galaxy has K a (total) = 12.86 mag in the 
2MASS Extended Objects catalogue. 
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(a) PKS 0022-60 (b) PMNJ 2008-6110 



(c) MRC 2041-617 (d) ESO 075-41 



(e) PKS 2356-61 (f) PKS 2356-61 with SUMSS 843 MHz contours overlaid 

Figure Al. Radio contours overlaid on SuperCOSMOS optical images of sources detected at 18 GHz and discussed in the text. 
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Figure A2. Radio contours overlaid on SuperCOSMOS optical images of sources detected at 18 GHz and discussed in the text. 
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(a) PKS 1105-680 (b) PKS 1133-681 



(c) PKS 1420-679 (d) PKS 1448-648 



(e)WKK5585 (f)NGC6328 

Figure A3. Radio contours overlaid on SuperCOSMOS optical images of sources detected at 18 GHz and discussed in the text. 
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Figure A4. Radio contours overlaid on SuperCOSMOS optical images of sources detected at 18 GHz and discussed in the text. 
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(a) PKS 2300-683 (b) PKS 2303-656 



(c) IRAS 23074-5957 

Figure A5. Radio contours overlaid on SupcrCOSMOS optical images of sources detected at 18 GHz and discussed in the text. 
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